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ABSTRACT 

The Great Observatories All-Sky LIRG Survey (GOALS) is a comprehensive, multiwavelength study 
of luminous infrared galaxies (LIRGs) in the local universe. Here we present low resolution Spitzer 
IRS spectra covering 5-38 /zm and provide a basic analysis of the mid-IR spectral properties observed 
for nearby LIRGs. In a companion paper, we discuss detailed fits to the spectra and compare the 
LIRGs to other classes of galaxies. The GOALS sample of 244 nuclei in 180 luminous (10'^"'^ < 
Lir/Lq < 10^^) and 22 ultraluminous {Ljji/Lq > 10^^) IR galaxies represents a complete subset 
of the IRAS Revised Bright Galaxy Sample and covers a range of merger stages, morphologies and 
spectral types. The majority (>60%) of the GOALS LIRGs have high 6.2 /xm PAH equivalent widths 
(EQW6.2/jm > 0.4 /xm) and low levels of silicate absorption (59.7^™ > -1-0). There is a general trend 
among the U/LIRGs for both silicate depth and mid- infrared (A4IR) slope to increase with increasing 
L/i?- U/LIRGs in the late to final stages of a merger also have, on average, steeper MIR slopes and 
higher levels of dust obscuration. Together, these trends suggest that as gas & dust is tunneled towards 
the center of a coalescing merger, the nuclei become more compact and more obscured. As a result, 
the dust temperature increases leading also to a steeper MIR slope. The sources that depart from 
these correlations have very low PAH equivalent width (EQW6.2/:im < 0.1 /xm) consistent with their 
emission being dominated by an AGN in the MIR. These extremely low PAH equivalent width sources 
separate into two distinct types: relatively unobscurcd sources with a very hot dust component (and 
thus very shallow MIR slopes) and heavily dust obscured nuclei with a steep temperature gradient. 
The most heavily dust obscured sources are also the most compact in their MIR emission, suggesting 
that the obscuring (cool) dust is associated with the outer regions of the starburst and not simply a 
measure of the dust along the line of sight through a large, dusty disk. A marked decline is seen for 
the fraction of high EQW (star formation dominated) sources as the merger progresses. The decline is 
accompanied by an increase in the fraction of composite sources while the fraction of sources where an 
AGN dominates the MIR emission remains low. When compared to the MIR spectra of submillimeter 
galaxies (SMGs) at z^2, both the average GOALS LIRG and ULIRG spectra are more absorbed at 
9.7 fim and the average GOALS LIRG has more PAH emission. However, when the AGN contributions 
to both the local GOALS LIRGs and the high-z SMGs are removed, the average local starbursting 
LIRG closely resembles the starburst-dominated SMGs. 
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1. INTRODUCTION 

A principal achievement of the Infrared Astronomical 
Satellite (IRAS) was the discovery of a large popula- 
tion of galaxies whose bolometric luminosities were dom- 
inated by emission in the infrared. At the highest lumi- 
nosities, local ultraluminous infrared galax ies (ULIRGs; 
L/_R > lO^^L©) have been heavily stu died ( lArmus et al 
20071 [Sanders et al.| |19 88 ; Murphy ct al. IMTT Spoon 
et al.|l20061 | Desai et al 2007. Rigopoulou et al.,,1999 



tra for 244 galaxy nuclei in the 202 nearby GOALS 
U/LIRG systems taken with the low resol ution module 



Genzel et al. 1998p , anoCa clear formation picture has 
been pieced together to explain their extreme emission 
in the infrared: more than 90% of local ULIRGs are the 
products of major mergers between molecular gas-rich 
galaxies. The large amounts of gas that are funneled 
into the centers of these mergers lead to intense star for- 
mation, the feeding of a central AGN, extremely compact 
reservoirs of molecular gas, and infrared luminosities on 
the order of ten times their optical luminosities. 

While ULIRGs constitute onl y 3% of the IRAS Rev ised 
Bright Galaxy Sample (RBGS; 
just slightly lower luminosities, 1 
ies (LIRGs; 10"M(t, < Ljr < lO^^Mr. 



Sanders et al. 2003), at 
uminous infrared galax- 



^ ^IR ^ iu'"ivi0j make up almost 
1/3 of the IR sources and have formation histories that 
are far less straightforward. In the local universe, there 
is evidence that galaxy-galaxy interactions and mergers 
drive the l arge IR luminosities in some LIRGs (Sanders 
& Mirabel 1996) and many high-z submillimeter galaxies 
(SMGs) s how hints of disturbed optical and rad io mor- 

. How- 
LIRGs 



phologies (Blain et al. 



ever, at least 20% and as many as 40% o: 



2002 Dasyra et al. 



2008)) 

ioca 



may have no history of major tidal interactions (How- 
ell et al., in prep). LIRGs are also represented across 
the full range of merger stages, unlike ULIRGs which are 
almost always at the very end stages of coalescing. 

Although LIRGs are relatively rare in the local uni- 
verse, their comoving number density increase s by more 



than 1 00 times fro m the cur rent epoch to z^l, (Le Floc'h 
et al.pOOS, MagneUi et al.|r2009) until LIRGs dominate 
the total IR energy density at redshifts of z~l-2 when 



star format ion in the universe was at its peak (Caputi 
et al.|[2007 ) . Piecing together the formation mechanisms 



and subsequent evolution of these LIRG systems is thus 
vital to understanding the processes governing star for- 
mation and black hole accretion, the main sources of 
emitting power in the IR. 

The Great Observatorie s All-sky LIRG Survey 
(GOALS; Armus et al.||2009 ) represents a complete sub- 
set of the RBGS comprising 180 LIRGs and 22 ULIRGs 
and aims to provide a multiwavelength understanding of 
the formation and evolution of local LIRGs as a class of 
galaxy. As part of the Spitzer Legacy survey, a complete 
set of IR imaging (Infrared Array Camera (IRAC) at 
3.6, 4.5, 5, and 8/im, and Multiband Imaging Photome- 
ter (MIPS) at 24, 70, and 160 /im) and IR spectroscopy 
at both high and low resolution (Infrared Spectrograph 
(IRS) from 5-38 /im) is available for the entire sample. In 
addition, imaging in the near- IR/optical (Hubb le Space 
Telescope NICMOS and ACS; |Haan et al.|[20TT| Kim et 
al., in prep) , the UV (Galaxy E volution Explorer near- 
and far-UV; [HoWell et al.|[2010l ), and the X-Ray (Chan- 
dra; Iwasawa et ai. 2UTT| rare available for large subsets 
of the sample. 
In this paper, we present the mid-infrared (MIR) spec- 



on th e Spitzer Infrared Spectrograph (IRS; [Houck et al 
20041. The MIR properties derived directly from such 



a large, complete sample of LIRG spectra will allow us 
to place these intermediate-luminosity systems into the 
context of both the extensive previous local ULIRG stud- 



starbursting systems (Brandl et al. 


2006 


Smith et al. 


2007b 


O'Dowd et al.||:^009 


Wu et al 


12010 


. 



and dust features as well as the SEDs covered by the IRS 
data, for the entire sample along with the comparison 
of MIR galaxy properties to those at other wavelengths 
will be presented in Stierwalt et al. (2013b). The analy- 
sis presented here focuses on properties derived directly 
from the MIR spectra. In Section [2) we present the low 
resolution MIR spectra observed with the Spitzer IRS 
Short-Low and Long-Low modules and describe our data 
reduction methods. In Section [3] we give the distribu- 
tions of the MIR properties and investigate correlations 
with L/fl and compactness. In Section |4) we follow each 
MIR property through the merging process, and we place 
our results into a high redshift context through compar- 
isons to MIR spectra of submillimeter galaxies at z~2 
in Section [5) Finally, our summary and conclusions are 
presented in Section [6) 

2. OBSERVATIONS & DATA REDUCTION 

2.1. The Sample 

The GOALS sample consists of 244 galaxy nuclei in 
180 luminous and 22 ultraluminous nearby IR galaxies. 
New spectra were obtained using the staring mode for the 
IRS Short-Low (SL: 5.5-14.5 fim) and Long-Low (LL: 14- 
38 /im) modules for 157 galaxies (PID 30323; PI L. Ar- 
mus). Integration times were determined from nuclear 
flux densities measured from IRAC and MIPS images 
and range from 45-120 seconds in SL and 30-120 sec- 
onds in LL. Secondary nuclei were targeted when the 
MIPS 24 fim flux ratio of primary to secondary was < 
5. Archival spectroscopic observations were used for the 
remaining 45 systems and borrowed most heavily from 
staring program PIDs 105, 3247, & 20549 and mapping 
program PIDs 73, 3269, & 30577. 

All 202 systems are nearby but cover a range of dis- 
tances (15 Mpc < D < 400 Mpc) and so the resulting 
projected IRS slit size varies from source to source. At 
the median galaxy distance of 100 Mpc, the nuclear spec- 
trum covers the central 1.8 kpc in SL and the central 5.2 
kpc in LL. 

2.2. Data Reduction 

Staring mode spectroscopic data were reduced using 
the S17 and S18.7 IRS pipelines from the Spitzer Science 
Centeipj For most sources, off-source nods were used 
to perform background sky subtraction. In the cases of 
more extended objects, dedicated background pointings 
were used to determine the sky surface brightness. One 
dimensional spectra were extracted using the standard 
extraction aperture and point source calibration modes 
in SPICErMwhich employs a tapered extraction aperture 



http:/ /ssc. spitzer. caltech.edu/irs/features. html 
http:/ /ssc. spitzer. caltech.edu/postbcd/doc/spice.pdf 
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that averages roughly to a size of 10"6 x 36'.'6 in LL and 
3'.'7 X 9'.'5 in SL. After masking bad pixels, multiple nods 
were averaged to produce the final spectrum. 

Of the archival data, 27 spectra were taken in staring 
mode and were reduced as described above. For the re- 
maining 18 systems, spectra were extracted f rom low res- 
olution mapping mode data using CUBISM ( Smith et al. 
2007a). Two-dimensional BCDs were asseinbled, obvi- 
ous bad pixels were removed and nuclear spectra were 
extracted. In two cases (CGCGOll-076 and IC1623B), 
smaller apertures were necessary to avoid other sources 
in the Long-Low maps, but for most sources 2x5 pixel 
extraction apertures centered on the galaxy's nucleus 
were used to resemble as closely as possible the results 
that would have been achieved with staring mode ob- 
servations. However, since the tapered aperture used by 
SPICE cannot be completely reproduced by the square 
apertures in CUBISM, a further mapping-to-staring- 
mode correction was applied to all spectra derived from 
low resolution archival maps. The correction, a mul- 
tiplicative factor that is a function of wavelength, was 
derived from NGC6240, a star-forming merger remnant 
typical of the GOALS sample for which both staring and 
mapping data were obtained. The correction function 
varies from 1.3 to 2.7 over SL wavelengths and from 1.7 
to 2.3 over LL wavelengths. 

The IRAC 8 /xm (Channel 4) images for six exam- 
ple GOALS galaxy nuclei are shown in Figure [T] with 
the SL and LL extraction aperture projections overlaid 
(in the case of staring mode data) or with the CU- 
BISM extraction windows overlaid (in the case of map- 
ping mode data). The low resolution IRS spectrum for 
each source is also presented along with each MIR image. 
Spectra for the remainder of the galaxy nuclei, ordered 
by right ascension, are available as online material and 
can also be found at http://goals.ipac.caltech.edu/. For 
five galaxies (IIIZw035, IRASF03359-I-1523, MCG-l-08- 
18-013, IRASF17132-h5313, and MCG-01-60-Q22), the 
archival SL staring mode observations were not cen- 
tered on the galaxy nucleus, so the SL slit overlays 
are not shown and the extracted spectra were not used 
in our analysis. Complete IRS observations were not 
available for an additional 6 galaxies (no LL spec- 
tra: NGC2388, NGC4922, and VV705; no SL spectra: 
IRASF08339+6517; no IRS dat a: ESO550-IG025 and 
IC4518). One galaxy (NGC1068; [Howell et al.|2007| ) sat- 
urates the spectrograph and so is also not shown. 

2.3. Scale Factors 

For each spectrum a break occurs between the SL and 
LL modules near 14 /im due to the larger LL slit, which 
covers nine times the area covered by the SL slit. The 
scale factors required to match the SL flux to the LL flux 
are not applied to the spectra in Figures[T]and Al but are 
calculated from the overlap in the SLl and LL2 modules 
and presented in Table [T] Scale factors are not given for 
any source missing either SL or LL data. For a small mi- 
nority of cases (CGCG448-020, ESO077-IG014, ES0173- 
G015, ESO255-IG007, ESO343-IG013, ESO440-IG058 
(northern nuclei only), IRAS03582-I-6012, IRASF06076- 
2139, NGC5653, NGC6090, NGC3690 (western nuclei 
only), and NGC5256), the scale factor is not given be- 
cause the placement of the LL slit covered multiple nuclei 
while the smaller SL slit covered only one. 



The median scale factors are 1.22 and 1.70 for the star- 
ing and mapping mode data respectively. The larger me- 
dian scale factor for mapping data most likely reflects a 
selection bias toward mapping more extended sources. 
Twelve scale factors are < 1 (i.e. more flux is recov- 
ered from SL than from LL), but for all twelve, the scale 
factors are also >0.9 and thus represent normal statisti- 
cal scatter for sources with scale factors near unity. No 
clear correlation is observed between the scale factors 
and galaxy distance, but at distances > 300 Mpc, a cut- 
off that includes 6 sources, the scale factors are all < 
1.2. Similarly, at distances closer than 30 Mpc, there are 
three GOALS sources that all have scale factors > 1.6. 

The scale factors are applied as a uniform multiplica- 
tive factor across the entirety of the SL spectra and thus 
boost equally the PAH fluxes, the continuum, and the 
absorption features. Since calculations of the equivalent 
width of the 6.2 /im PAH and the depth of the silicate 
feature at 9.7 fim. (EQW6.2;i?n and sg,r^m', see next sec- 
tion) both use measurements of feature flux relative to 
the continuum, neither are affected by the scaling of the 
SL spectrum at these low redshifts. The MIR slopes 
(Fi,[30/im]/F,y[15/im]) are also unaffected as they only 
rely on data from the (unsealed) LL portion of the spec- 
trum. 

2.4. S9.7^m, MIR Slope, & EQWe .2/im 

Silicate depths at 9.7 /im (sg.y^^m) were measured di- 
rectly from the MIR spectra via: s\ = log{f\/C\) where 
fx is the measured flux at the central wavelength of the 
absorption feature and C\ is where the level of the con- 
tinuum flux would be in the absence of the absorption 
feature, based on an extrapolation to the surrounding 
continuum. Thus, a positive value, s\ > 0, suggests 
emission at that wavelength and the deeper the absorp- 
tion, the lower the S9.7^m value. 

The fluxes F^ at 15 /zm and at 30 /xm were determined 
from the average of eight data points surrounding each 
wavelength and were then used to calculate the MIR 
slope. The wavelength regions used fell within ^14.7- 
15.4 /im for F^[15fim] and -29.5-30.8 /im for F^[30/im]. 

Equivalent widths for the 6.2 /im PAH feature 
(EQW6.2/im) were measured for each sp ectrum using the 
method outlined in Brandl et al. (2006). Briefly, a spline 
fit was used to estimate the continuum surrounding the 
6.2 /im PAH feature, and the continuum fit was sub- 
tracted from the spectrum. In most cases, anchor points 
in determining the continuum were set at 5.15/xm < 
A <5.31 /im, 5.8 /im < A < 5.9 /im, 6.5 /xm < A < 6.8 /xm, 
and 7.1 fim < A < 7.2 /tm, but each spectrum was visually 
inspected to make sure no features or bad points occurred 
in these ranges. The PAH flux was then measured using 
direct integration. The 6.2 /im feature was selected for 
the EQW calculation because, of the five brightest PAH 
features, it is the least affected by silicate absorption at 
9.7 /im and 18.5 /im, and it is not blended with other PAH 
features. However, in some cases the 6.2 /im PAH feature 
partially overlaps with the absorption feature due to wa- 
ter ice at 6.0 /im. For those sources found by the spec- 
tral fitting to have r^ce >0 (see Stierwalt et al. 2013b), 
the ice absorption was assumed to affect the underly- 
ing continuum but not the PAH emission, and the EQW 
was calculated accordingly. Four galaxies have only up- 
per limits placed on their EQW6.2,.m: IRAS05223-M908, 
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MCG-03-34-064, NGC4418, and IRAS08572+3015. 

2.5. Merger Stages 

Merger stages for the entire sample were determined 
via visual inspection of the IRAC 3.6 /Ltm (Channel 1) im- 
ages. Each galaxy was assigned one of the following five 
designations: 'N' for nonmergers (no sign of merger ac- 
tivity or massive neighbors), 'a' for pre-mergers (galaxy 
pairs prior to a first encounter), 'b' for early-stage merg- 
ers (post-first encounter with galaxy disks still symmetric 
and in tact but with signs of tidal tails), 'c' for mid- 
stage mergers (showing amorphous disks, tidal tails, and 
other signs of merger activity), or 'd' for late-stage merg- 
ers (two nuclei in a common envelope). Given the res- 
olution of the IRAC images (~2"), late stage mergers 
can be easily mistaken for nonmergers in the 3.6- /xm im- 
ages. To alleviate this problem, any galaxies classified as 
nonmergers or early stage mergers in the IRAC images 
with available higher resolution imaging in the literature 
that clearly showed signs of a late stage major merger 
were changed accordingly. We also use the literature to 
identify spectroscopic pairs which resulted in reclassify- 
ing some nonmergers as pre-mergers. 

For a subset of 78 GOALS galaxies (all with 
log(L/fl;/LQ) >11.4), we have additional merger classifi- 
cations based on available HST B, I, and H-band images. 
The higher resolution of this imaging enables a more de- 
tailed classification system with more finely tuned merger 
stage designations (stages through 6). These merger 
stages were already described and presented in |Haan| 
et al. (2011 1, but we reproduce and discuss them here to 
aid with cross-referencing the two classification schemes. 
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3. MID-INFRARED PROPERTIES OF NEARBY LIRGS 

3.1. LIRG vs ULIRG Distributions 

Silicate depths, MIR slopes, PAH equivalent widths, 
and all associated uncertainties for the GOALS sample, 
in addition to the SL-to-LL scale factors and merger 
stages, are presented in Table [l] and the distributions 
of EQW6.2/^m,S9.7^m, and MIR slope are shown in Fig- 
ure [2] The EQW6.2/im and Sg y^^ parameters are not 
given for the five sources with off-centered SL spectra, 
and MIR slopes are not presented for the four sources 
without available either SL or LL spectra or for the 12 
sources for which multiple nuclei are observed within the 
LL slit. 

As shown in Figure [21 the majority of LIRGs (63%) 
are dominated by PAITemission (EQW6.2/^m >0.4/im), 
show little to no silicate absorption (sg.y^m > -1), 
and have MIR slopes of 4 < F^ [30/.tm]/Fi,[15 ^m] < 
10. Only six LIRGs have deep silicate absorption 
with S9 7^m < -1.75 (NGC4418, IRAS03582-f6012_E, 
ESO203-IG001, IRASF10038-3338, IRASF12224-0624, 
and ESO60-IG016). The remainder of the LIRGs show 
weak to no silicate absorption with a significant fraction 
(23%) of LIRGs showing silicates in emission at 9.7 /im, 
including 11% with sg y^m > 0.15. A few of the LIRGs 
with S9.7^„i > are likely AGN-dominated (EQW6.2pm 
< 0.27 yum) and thus any absorption at 9.7 /xm is filled 
in by an excess of hot dust. However, most are lower 
luminosity galaxies with 90% having log(L/fl/L0) < 
11.25. These LIRGs are likely analagous to the unob- 
scured starburst NGC 7714, a galaxy whose IR emis- 
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Fig. 2. — Distributions of MIR spectral parameters (upper pan- 
els). Top: silicate strength at 9.7/im, Middle: logarithm of MIR 
slope, and Bottom: equivalent width of the 6.2 /xm PAH feature. 
On average, GOALS ULIRGs (filled red histograms) have deeper 
silicate absorption depths, steeper MIR slopes, andlower equivalent 
widths than the GOALS sample as a whole (white histograms). 
The lower panel on each plot shows the same GOALS ULIRG dis- 
tributions with a smaller y-scale. The filled solid black portion 
of the lowest bin of the EQW6.2(jm histogram represents the four 
sources for which only upper limits are measured. 



sion is fueled alm ost entirely by star formation (Mar- 
shall et al. 2007). The silicate strengths in the LIRGs 
have a median of sg.r^m = -0.25 ± 0.58 and range 
from the heavily obscured NGC4418 at sg y^^ = -3.51 
± 0.09 to NGC5395, the southern component of the 
LIRG system Arp84, which shows silicates in emission 
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(sg.y^m = 0.52 ± 0.07). Five LIRGs are continuum dom- 
inated and show at most only weak PAH or line fea- 
tures (EQW6.2,™ < 0.04 ^m and 59.7^™ > -0.2; MCG- 
03-34-064, IRAS05223-t-1908, NGC1275, NGC7674, and 
AM0702-601_N). 

While the majority of LIRGs favor the high end 
of the distribution in both EQW6.2;im and sg.T^im, 
they are found clustered in an intermediate range of 
MIR slopes with a median of Fi,[30/im]/Fy[15^m] — 
7.11±4.74. The MIR slopes measured for the LIRGs 
range from F^,[30Aim]/F^[15/im] = 2.00 ± 0.01 in 
IRAS05223-I-1908 which shows a near power-law spec- 
trum in the MIR to F^[30Aim]/F^[15Atm] = 35.40 ± 1.38 
in IRAS10173-t-0828. 

For those LIRGs with measurable 6.2 ^m PAH EQWs, 
the values range from EQW6.2^im — 0.005 /im ± 
0.003 /im for the northeastern component of the LIRG 
pair IRAS03582-H6012 to EQW6.2,™ = 0.78 /im ± 
0.01 /im for the most southeastern of the three galax- 
ies composing the LIRG system IRAS17578-0400. The 
distribution for all of the GOALS LIRGs has a median 
of EQW6.2pm = 0.55 /im ±0.18/im. The same median 
value was found f or a sample of lowe r luminosity star- 
bursting galaxies (Brandl et al. 20061. Tight limits are 
placed on the EQ W tor the three LlKGs and one ULIRG 
without a 6.2 /im PAH detection: IRAS05223-1-1908 at 
<0.043/im, MCG-03-34-064 at <0.044/im, NGC4418 at 
<0.066 /im, and IRAS08572+3915 at <0.081 /im. 

The GOALS ULIRGs, represented by the solid red his- 
tograms in Figure [2j show a clear offset from the LIRGs 
in their distributions for all three fundamental proper- 
ties. The ULIRGs have a higher median flux density ra- 
tio (F;,[30/im]/Fi.[15/im] = 12.54±5.41), a lower median 
PAH equivalent width (EQW6.2Mm = 0.30 /im±0.17 /im), 
and deeper median silicate absorption (sg — -1.05 ± 
0.85). The GOALS ULIRGs span nearly the full range 
of MIR slopes covered by LIRGs but are not found with 
EQW6.2pm >0.52/im or with S9.7;^m >-0.15. Comparing 
the derived values fo r the 22 ULIRGs in G OALS with the 
larg er samples from ISpb on et al.| (104 ULIRGs) 

and IVeilleux etal] ([2TO) CQ^UEST; 50 ULIRGs), we find 
that the larger numbers of ULIRGs in these samples re- 
sult in a larger spread in MIR properties (i.e. 6.2 /im 
PAH EQWs up to .8 /im and silicate depths up to 0.2; 
Spoon et al. (2007)). However, the median values are 
consistent with ULIRGs having lower EQW6.2^im, deeper 
sg.7fj,rm and steeper MIR slope than LI RGs: median 
EQW fi 2i/,iT7, = 0.15 /im & Sg 7^„i = -1.47 dSpo on et al 
20071) and median Fi.[30/im]/Fi,[15/im] 
efaL][2009l ). 
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I'he results of a Kolmogorov-Smirnov (KS) test give 
probabilities of <0.01% that the chance deviations be- 
tween the distributions of EQW6.2/im, sg.T/im, and MIR 
slope for GOALS LIRGs vs ULIRGs are expected to be 
larger assuming they are derived from the same parent 
sample. In other words, the two samples are significantly 
different. These probabilities decrea se by several orders 
of magnitude when the QUEST and Spoon et al. (20071 
ULIRGs are included. When the GOALS ULlRGsl^re 
compared to the larger ULIRG samples, the KS test 
suggests the chance deviations in their distributions in 
MIR slope, EQW6.2Aim, and sg.7^m are expected to be 
larger with probabilities of 80%, 4 0%, and 30% , i.e. it is 
likely the GOALS ULIRGs and thelSpoon et al.|(l2007l & 



Veilleux et al. ( 2009 ) samples are derived from the same 
parent sample. 

3.2. Correlations with Lm 

Figure [3] shows the distributions of sg.7^m, MIR slope, 
and EQW6.2/im as a function of IR luminosity, Ijir. The 
IR lumino sities for all 202 U/LIRG systems were pre- 
sented in Armus et al.' (2009) and derived using the defi- 
nitions of [Sanders fc Mirabel ( ,1996p | In cases of multi- 
ple nuclei, the total L/^j tor tne system is divided accord- 
ing to the ratio of the fluxes at 70 fim for each nuclei. In 
a small number of cases, 70 /im images are not available 
and so 24 /tm flux ratios are used instead. 

There is a general trend among the U/LIRGs for both 
silicate depth and MIR slope to increase with increas- 
ing L/fl. The sources that depart from these correla- 
tions at deep levels of silicate obscuration (top panel) or 
shallow MIR slopes (middle panel) have, in both cases, 
very low PAH equivalent width (EQW6.2/j?ri < 0.27 /im) 
and are thus likely dominated by emission from an AGN. 
Increasingly luminous systems become increasingly dust 
obscured until a turnover occurs at 59.7^™^ -1.5, above 
which the buried AGN candidates show no further cor- 
relation between sg.7pm and L//j. As hjR decreases, the 
MIR slopes flatten until F^[30/im]/Fi.[15/im]< 0.5, be- 
low which the relatively unobscured AGN have high Lm 
given their slopes. ULIRGs have an average EQW6.2/im 
that is lower than that for LIRGs, but sources with a 
large range of luminosities are found at each equivalent 
width (lower panel) so there is not a tight correlation 
between hjR and EQW6.2;^m- 

3.3. Disentangling sg,Tfj,m, MIR Slope, & EQWe,2fj,m 

To further disentangle the relationship between the 
3 main MIR parameters, we examine the Sg.7^m and 
EQW6.2/im versus MIR slope parameter spaces in Figure 
[4[ The distribution of sg.7^TO with MIR slope is color- 
coded by EQW6.2/im (panel a) while the distribution 
of EQW6.2Atm with MIR slope is color-coded by sg.T^rn 
(panel h). 

In the majority of LIRGs, star formation dominates 
the intense MIR emission and the similar conditions in 
the photodissociation regions resulting in this (PAH- 
dominated) emission lead to similar MIR properties 
among the bulk of the GOALS sample. In Figure |4^, 
the majority of GOALS galaxies, those with EQW6.2/im 
> 0.27 /im (green squares and blue stars), show a rough 
correlation between increasing MIR slope and increasing 
silicate depth (lowe r sg.7^m). The starburst galaxies of 
Brandl et al] ( |2006[ ) which span mostly luminosities be- 
low IO^^Lq exhibit a similar relationship between rg 7^^ 
and Fi,[30/im]/Fiy[15/im] with the same slope. The trend 
in Figure suggests that the average dust temperature 
rises as a consequence of the nuclei becoming more ob- 
scured and compact. As the dust temperature increases, 
the rising portion of the blackbody emission spectrum 
shifts to shorter wavelengths, and warmer sources have 
increasingly more flux at 30 /xm as seen for the GOALS 
U/LIRGs with sg.7^m >-1.5. 

22 Lir/Lq = 47r(Di,[m])2(F/fl[Vym-2])/3.826 X 10'^<^[Wm-^] 

and FfR = 1.8 X 10-"(13.48/i2^m + 5.16/25^™ + 2.58/60^™ + 
/i00Mm[VKm-2]) 
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Fig. 3. — Distribution of MIR spectral parameters with hm 
color-coded by EQW6.2pm. Top: siUcate strength at 9.7/xm, Mid- 
dle: logarithm of MIR slope, and Bottom: equivalent width of 
the 6.2 fj,m PAH feature. There is a loose trend among LIRGs for 
increasing silicate depth and MIR, slope with increasing hjn. How- 
ever, LIRGs span nearly the full range of EQW6.2fim at any given 
luminosity. 



Most of the sources with low PAH equivalent width, 
however, do not follow these simple trends in MIR prop- 
erties. In Figure l4k, these low-EQW sources (red cir- 
cles) are split roughly into two populations: those that 
are relatively unobscured with shallow MIR slopes and 



those heavily obscured sources (s9.7;im<-l-5) with steep 
MIR slopes. A similar split is observed in Figure [4]d: 
for EQW6.2/im <0.27^m, the heavily obscured sources 
(purple circles) are found at steeper MIR slopes while 
the relatively unobscured sources (magenta circles) are 
found at the shallowest flux density ratios. 

An increasingly significant hot dust component from 
an AGN leads both to a decrease in EQW6.2;im and 
to a flatter MIR slope. For the GOALS sources with 
the strongest, relatively unobscured AGN (s9.7^m>-0.6; 
EQW6.2;im^0.05 /im), an upper limit to the MIR slope 
can be set at F,y[30/im]/Fj,[15/im]<4 from both panels 
in Figured These galaxies (including IRAS05223-I-1908 
and UGCD8058) are represented by the red circles in the 
lower right corner of panel a and the magenta circles in 
the lower left corner of panel b. No other sources are 
found with flatter MIR slopes. This li mit agrees with 



that observed for the starbursts of Brandl et al. ( 2006 1 



and for the QUEST ULIRGs of Veilleux eTaTrpnogp ^. 
Relatively unobscured AGN can thus be identihed based 
on their low MIR flux density ratio alone. 

In the most heavily obscured, low EQW galaxies, how- 
ever, the MIR continuum slopes are steeper due to the 
buried, hot source. These galaxies (red circles in the 
left half of Figure |4^ and purple circles in Figure |4jD) 
have steep MIR slopes for the same reason sources with 
S9.7;im~-1.5 in Figure |4^ have steep MIR slopes: most of 
the warm dust emission is hidden behind a large amount 
of cooler dust. A comparison to the 5 mJy Unbiased 
Spitzer Extragalactic Survey (5MUSES; [Wu et al.||2010 | 
highhghts the difference between the low and high sg.r^^m 
sources at low EQW6.2Aim- The 5MUSES sample is 24- 
fim selected (indicating the presence of hot dust) but 
lacks the heavily obscured sources found in GOALS. The 
distributions for the two samples in FigureHjD are roughly 
the same (5MUSES is represented by the dashed lin e al- 
though ther e is significant scatter about this line; see |Wu| 
et al. (2010)) - both show the lowest EQW6.2Aim sources 
have the shallowest MIR slopes - except 5MUSES lacks 
the obs£ured low EQW6.2/^m galaxies (purple circles in 



Figure |4|: 

The apparent strength of the 9.7 //m silicate feature 
(i.e. the depth of the absorption feature that does not 
account for any silicate emission, sg.y^^m) is shown versus 
EQW6.2Aim in Figure js] for GOALS LIRGs (open circles) 
and ULIRGs (red triangles). No galaxies are observed 
with both high equivalent widths and large levels of sili- 
cate absorption. However, at low equivalent widths, two 
distinc t branches, similar to those seen by [Spoon et al.| 
( 2007 1 , emerge that clearly distinguish the lower equiva- 
lent width (EQWe,2^m < 0.1 /zm) sources with minimal 
to no silicate absorption (sg.r^^m > -0.5) from those dom- 
inated by silicate absorption (heavily obscured sources; 
sg.7fj.m <-1.75). Sources with intermediate levels of sili- 
cate absorption are not found at low equivalent widths. 

As shown in Figure [Sj the highly absorbed sources are 
not limited to ULIRGs. At values of sg.y^m < -1.75, the 
GOALS sample includes five ULIRGs (labeled in Fig- 
ure p| as well as the dense, compact nascent starburst 
LIRU NGC4418 (ISpoon et al. 2001; Roussel et al.||20Q3" 



Evans et al.| 2006 1 and five additional LIRGs thai; span 
a large range of LIRG luminosities: IRASF12224-0624 
(log(L/7?,/L0) = 11.36), IRAS03582+6012 (log(L//?/L0) 
= 11.42), IRASF10038-3338 (log(L/ii/L0) = 11.78), 
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Fig. 4. — Distribution of MIR slope (F^ [30 /xm]/Fi, [15 /im]) versus a) silicate absorption at 9.7 /^m color-coded by EQW6.2/jm and b) 
EQW6.2pm color-coded by sg.T^m. GOALS sources with EQW6.2(/m >0.27/xm (green squares + blue stars in panel a, right side of panel 
b) show a rough correlation between increasing silicate depth and increasing MIR slope (a) and follow the correlation between EQW6.2(jm 
and MIR slope observed in the 24- /xm selected 5MUSES sample (dashed line, b; Wu ct al. (2010)). At low EQW (EQW6.2Mm <0.27Atm; 
red circles in panel a, left side of panel b), relatively unobscured AGN-dominated sources all have MIR slopes below ~4. However, at 
the deepest levels of silicate absorption (left side of panel a, purple circles in panel b), the MIR slope is no longer a clear indicator of 
temperature and so the heavily obscured sources do not follow the trend in panel a, and the location of the 15- /xm continuum between the 
9.7/im and 18.5 /^m absorption features lead to elevated MIR slopes in panel b. Although far less numerous (only 18% of GOALS nuclei 
have EQW6.2;tm < 0.27/jm), the lowest equivalent width sources cover a wider range of L/ij, MIR slope, and S9.7(jm than those sources of 
higher EQW6.2fjm that make up the majority of the sample. 
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Fig. 5. — Equivalent width of the 6.2 ^m PAH versus silicate strength for the LIRGs (open circles) and ULIRGs (red triangles) of the 
GOALS sample. The majority (>60%) of the LIRGs are found at EQW6.2nm > 0.4/im and S9.7^m > -1.0, while nearly all ULIRGs 
have EQW6.2nm < 0.5/^m and sg.T^m < -0.5. Sources at low EQW are divided into two distinct branches (silicate-dominated versus 
continuum-dominated) . 
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ESO60-IG016 (log(L/7i/Lo) 
IGOOl (log(L/^/L0) = 11.86). 



3.4. 



11.82), and ESO203- 



Compactness 

The most heavily obscured nuclei among the GOALS 
galaxies are also the most compact: they all show little- 
to-no MIR emission extending outside of the IRS slit. 
In Figure |6j the silicate strength is plotted against rj, a 
parameter that represents the fraction of the emission at 
24 /im intercepted by the IRS slit: 



Zog(F*i^^^ [24^.™] /i^/;f/ [24Mm] ) 



(1) 



where F^{^^^ [24^m] is the total flux of a source as mea- 
sured from its MIPS 24 /xm image (Mazzarella et al., 
in prep) and F^iff [24/xto] is the flux within the IRS slit 
derived by convolving the MIPS-24/im flltcr with the 
low resolution IRS spectrum. The most obscured sources 
(sg.T/im < -1-75) all have 77 ^ (i.e. all of the flux mea- 
sured by the larger MIPS field of view at 24 /im is also 
recovered within the much smaller IRS slit). 

Although distance effects could act to disguise an ex- 
tended component in comparisons of total versus intra- 
slit fluxes if all of the obscured sources were the most dis- 
tant, the median distance for the heavily obscured, low 
r] nuclei is only 60% larger than the median distance for 
the less obscured sources (190 Mpc vs 115 Mpc) suggest- 
ing distance alone cannot be driving the difference in rj. 
Even more importantly, the heavily obscured, low ry tail 
of the distribution in Figure [6] includes the two closest, 
obscured LIRGs, NGC4418 at 36.5 Mpc and sg =- 
3.51 ± 0.09 and NGC3690 at 50.7 Mpc and sg.j^m =-1.65 
± 0.02. Both of these LIRGs would have been easily 
resolved had they shown any extended MIR emission. 
Additionally, the existence of galaxies with high 77 and 
low Sg.y^m with distances well above 100 Mpc indicate 
that extended sources can still be resolved even at larger 
distances. 

The fraction of resolved emission within the IRS slit is 
also much lower for the obscured sources. The fraction of 



extend ed emission (FEEi3.2;jjn) is deflned by Di'az-Santos 
et al.| (2010) as the fraction of emission withm the iRb 
slit origmatmg outside of the unresolved component (i.e. 
a point source at that distance). For the most obscured 
GOALS nuclei, the average (FEEi3.2/im) = 0.07 com- 
pared to (FEEi3.2;^m) — 0.39 for the remaining (weakly 
ob scured or unobscured) LI RGs. As discussed in detail 



Di'az-Santos et al. 



(2010), such a dramatic difference 
in fKhj between obscured and unobscured nuclei cannot 
be the result of distance effects alone. 

Together the low ry, the low FEEi3.2^m, and their in- 
clusion of nearby sources suggest the nuclei in these heav- 
ily obscured sources dominate the 24- /im emission from 
their parent galaxies, and so the most heavily obscured 
LIRGs and ULIRGs also have the most compact MIR 
continuum emission. Given their low EQW6.2/jm, if ex- 
treme levels of obscuration are not simply masking the 
PAH emission, the higher densities in these nuclei may 
create an environment where PAH dust grains are not 
present or the conditions are not appropriate for exciting 
them (i.e. lacking in photodissociation regions). Alter- 
natively, the low EQW6.2^m may indicate an increase in 
the continuum flux at 6 /im due to a rise in dust temper- 
ature. None of the low r/, high sg.r^m nuclei are observed 



to be [NeV] emitters ( Petric et al.pOll ), most likely be- 
cause their large optical depths obscure any line emission 
that would be present at 14.3 /^m. 



^ 1.5 



(N 



[I; 
(N 







0.5 



1 1 1 ' ' 

; .EQW,,^„,<0.27^m 


1 


0.27nm<EQW6,^„,<0.54nm 




<rEQW,,^„,>0.54jim 


-•- 




• 

ffl 




* 






1 1 1 < < 


9' 



-4 -3 -2 -1 
Silicate Strength i&^j^^J 



Fig. 6. — Silicate strength versus r], a measure of the total-to-sht 
flux ratio at 24 fim. GOALS LIRGs and ULIRGs are color coded 
by 6.2 lira PAH equivalent width with low equivalent width (AGN- 
dominated) sources (EQW6.2fjm < 0.27 fim) represented by red 
circles and high equivalent width (starburst-dominated) sources 
(EQWe. 2fim > 0.54 /^m) represented by blue stars. Intermedi- 
ate EQW (composite) sources are shown by orange (0.27 ^im < 
EQW6.2Mm < 0.41 /xm) and green (0.41 /xm < EQW6.2A.m < 
0.54 fj,m) squares. Heavily obscures sources have no extended com- 
ponent to their 24 /im emission {rj ~0). 



4. TRACING MIR PROPERTIES THROUGH MERGER 
STAGE 

In Figure [7] silicate strength, MIR slope, and PAH 
equivalent width are traced through merger stage for 
GOALS galaxies. To look for subtle differences in MIR 
properties throughout the merging process, we focus on 
only those sources that have HST classifications (column 
11 in Table [1]). Since this subset contains only six galax- 
ies with no mdication of a merger (stage 0), we include 
all of the nonmergers (Stage N) from the IRAC-based 
classifications (column 10 in Table [T]) to derive a more 
secure median for each spectral property. Although the 
HST data samples only LIRGs with log(L/fl/Lo)>11.4, 
the dense sampling of merger stages made possible by 
the deep, high spatial resolution optical and NIR images 
provides a much finer look at the spectral changes along 
the merger sequence. Mean values for each merger stage 
clipped at 3cr are shown in red with their associated stan- 
dard deviations. 

As the mergers progress and gas & dust is funneled 
towards the center, galaxies become on average more 
obscured with steeper MIR slopes. Silicate depths of 
sg.Tfim^ -1 are only reached at merger stages of 3 
and later. No LIRG systems in merger stage 1 have 
Fiy[30//m]/F,y[15/im]> 1, while the average MIR slope is 
>1 for the later stages 4-6. These two results agree with 
several studies finding higher L/^j at later merger stages 
since, as shown in Figure [3] increasing MIR slope and 
silicate depth are also linked to higher Lju in LIRGs. 
As merging galaxies coalesce, the nuclei become more 



10 



Stierwalt et al. 



compact and more obscured, and, as a result, the dust 
temperature increases leading to a steeper MIR slope as 
discussed in Section |3l 
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Fig. 7. — MIR properties of GOALS galaxies traced through 
merger stage. Top: silicate strength (sg.xy^m ), Middle: MIR slope 
(log(Fy[30/xm]/F^[15/im]), and Bottom: EQW6.2/jm- Mergers 
(Stages 1-6) are represented by the 78 GOALS galaxies for which 
high resolution HST imaging is available (Haan et al. 201 1[ ; see 
Column (11) in TablejlJ. Nonmergers (Stage JV) are classihed using 
IRAC 3.6 fim images and the literature (see Section [2.5| for details). 
Mean values for each merger stage clipped at 3ct are shown in red 
with their associated standard deviations. 



There is some indication that lower PAH equiva- 
lent widths are favored at later merger stages but this 
is mostly dominated by the fact that only starburst- 
dominated galaxies (EQW6.2;im>0.54 ^m) are observed 
in stage 1. (The one exception is southern component 
of the LIRG system AM0702-601.) For all other merger 
stages, the full range of EQW6.2/jm is observed. 



A clearer link between PAH equivalent width and 
merger stage is observed when galaxies are binned by 
their EQW6.2/^m(and thus the likely AGN contribution 
to the MIR). In Figure [s] the LIRGs are divided into 
three EQW6.2/ijn bins indicating AGN dominated sources 
(red circles), composite sources (green squares), and star- 
bursts (blue stars). Starbursts clearly play a domi nant 
role at ea r ly me rger stages as was also s hown by |Pet- 
ric et al.j ( |2011l ) and |Haan et all ( |2011[ ), but the de^ 
cline in the star burst contribution is not balanced by an 
increase in AGN-dominated sources. The contribution 
from LIRGs with an AGN dominating in the MIR stays 
at a roughly constant fraction throughout the merger 
process, but composite sources (i.e. the weaker AGN 
that are not yet entirely dominant over star formation in 
the MIR) show a marked increase at later merger stages. 
This may indicate that the timescales for the AGN to 
begin to dominate the MIR emission are longer than the 
merger timescale (a few hundred million years). In both 
Figures [7] and [Sj the nonmerging LIRGs cover nearly the 
full range of every MIR property investigated. 
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Fig. 8. — Overall AGN fraction (as determined by EQW of the 
6.2/im PAH) traced through merger stage for the GOALS sample. 
Merger stages are classified as described in Figure [7] A marked 
decline is seen for the fraction of high EQW (star formation domi- 
nated; blue stars) sources as the merger progresses. This decline is 
accompanied by an increased contribution not from the strongest 
AGN (red circles) which remain low but from the composite sources 
(i.e. weaker AGN that are not yet entirely dominant over star for- 
mation in the MIR; green squares). 



5. COMPARISONS TO SUBMILLIMETER GALAXIES 

The dust-enshrouded, strongly starbursting nature of 
LIRGs makes them obvious candidates for possible lo- 
cal analogs to the dusty submillimeter galaxies (SMGs) 
that make a significant contribution to the global star 
formation rate density at higher redshifts. In Figure [9j 
we compare different subsets of the GOALS MIR spectra 
with the ave rage SMG spectra from |Menendez-Delmestre] 
et al.| ( |2009[ ) (hereafter M09) derived from a sample of 24 
SMGs at redshifts of 0.65 < z <3.2. All average spec- 
tra for both the LIRGs and the SMGs are normalized at 
6. 



/irn . In agreement with the conclusions of |Desai et al.| 
(|2007| and M0_9, the average local ULIRG spectrum (red 
is more absorbed than the average SMG 
line) but has weaker PAH emission. Al- 



Ime m Figure K3r 
spectrum (black 
though GOALS LIRGs are less obscured than ULIRGs 
on average, the average LIRG spectrum (dashed line) 
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is still more absorbed than the average SMG while also 
showing stronger PAH emission at 6.2 /im, 7.7 ^m, and 
11.3 /xm. Even when the nuclear emission of galaxies is 
removed, the spectrum of the extended component of 
LIRGs does not resemble that of the total SMG com- 
posite ( Diaz-Santos et a"L||2011[). T he average local star- 
burst (blue line; Brandl et al.|2(J(J6 1 shows a similar level 
of silicate absorption but much stronger PAH emission 
compared to the average SMG. 




- b)AGN LIRGs 

- Starburst LIRGs — 
Starburst SMG (M09) 




Fig. 9. — Comparison of GOALS average LIRG and ULIRG 
spectra to avera ge submillimeter galaxy spectra from Mcncndcz- 
[Delmestre et aL] |2009)) : a) the composite SMG spectrum (blaclt) 
IS less obscured than the average ULIRG (red) and the average 
LIRG (dashed) but has weaker PAH emission than local starbursts 
(blue), b) after removing the AGN-dominated systems from both 
the average SMG and the average LIRG, the average starburst 
SMG spectrum (black) is well-represented by the starburst LIRGs 
(EQWe 2iim >0.54/^m; blue) but not the AGN-dominated LIRGs 
(EQWe. 2/im <0.27/im; red) with the exception of a feature at 
~10.5/im. All average spectra are normalized at A = 6.8 ^m and 
the shaded gray area represents the l-tr standard deviation to the 
averaged SMG spectrum. Residuals are shown in the bottom pan- 
els. 



The fraction of MIR emission attributed to AGN over- 



all for the GOALS LIRGs is only 12% ( |Petric etaL|2011[ ), 
and M09 observed a contribution of <32% from AGIN to 
the total bolometric luminosity in SMGs. However, as 
seen in Figures [l] and Al, those sources dominated by 
AGN have MIR spectra that are vastly different from 



those with strong PAH emission. The low equivalent 
width sources (i.e. those with MIR emission that is most 
likely AGN-dominated) also show a larger scatter in their 
MIR properties, as discussed in Section [3) To reduce pos- 
sible confusion caused by this AGN contribution, we also 
compare the average SMG spectrum for only those SMGs 
without AGN indicators in the MIR (i.e starburst SMGs 
with EQWy.y^m >1 A*m & aMiR <0.5; M09) to average 
LIRG spectra with and without an AGN contribution in 
Figure (9)3. 

The average AGN-dominated LIRG spectrum (LIRGs 
with EQW6.2pm <0.27/im; red line) clearly does not 
resemble the average starburst SMG (solid black line). 
However the average starburst-dominated LIRG spectra 
(EQW6.2/im >0.54 /.tm; blue line) is a better match to the 
average starburst SMG. All three average spectra agree 
within 15% below 10 /im (see the residuals in the lower 
panel) . 

Between 10-11 fjjn, all three average LIRG spectra in 
Figure [9]3 agree closely, but no subset of the GOALS 
LIRGs or ULIRGs reproduces the emission feature ob- 
served in the average starburst SMG spectrum near 
10.5 /im. Although the feature is also not detected in 
the average (low resolution) local starburst spectrum (see 
Figurejok), both the [SIV] emission line at 10.51 jjum and 
a PAH feature at 10.60 /im are clearly seen in the aver- 
age of the high resolution IR S spectra of the same star- 
burst sample (see Figure 4 oflBernard-Salas et al.|2009). 
The feature detected in the low resolutio n SMG Spec- 
tra is likely a blend of these two features (Sturm et al 
20001 iBernard-Salas et all 120091 ISmith et al.||2007bi' 



but may be dominated by the PAH feature emission 
since it remains faint and unresolved at low resolution. 



Bernard-Salas et al. ( 2009 1 also detect a third feature at 
10.75 /im that they associate with PAH emission due to 
its close correlation with the 11.3 /im PAH which may 
contribute to the emission in the SMGs. 

6. SUMMARY & CONCLUSIONS 

We presented low resolution IRS spectra for 244 galaxy 
nuclei in the GOALS sample of 180 LIRGs and 22 
ULIRGs. The GOALS galaxies cover a range of spec- 
tral types, silicate strengths, and merger stages, and rep- 
resent a complete subset of the IRAS Revised Bright 
Galaxy Sample. We investigated the MIR properties 
directly measured from the spectra and discovered the 
following: 

1) Local LIRGs cover a large range of MIR properties and 
any single LIRG cannot represent the class as a whole. 
LIRGs span 0.005 /im < EQW6.2Mm <0.78/im (with 
nondetections of the 6.2 /im PAH reaching 
EQW6.2Mm <0.043/im), -3.51 < < 0.052 
(with 23% of LIRGs showing silicate emission), and 2.00 

< Fj,[30/im]/F,^[15/i,m] < 35.40. However, the majority 
(63%) of LIRGs have EQW6.2A.™ > 0.4, sg.y^™ > -1.0, 
and MIR slopes in the range of 4 < Fy[30/im]/Fy[15/im] 

< 10. 

2) The GOALS ULIRGs span a narrower range of MIR 
properties than those covered by the LIRGs. When com- 
pared to LIRGs, the ULIRGs (L/^. > IO^^l©) have 
a steeper median slope (Fi,[30/im]/Fi,[15/im] = 12.54 
for the ULIRGs compared to Fi,[30/xm]/F,y[15/im] ~ 
7.11 for the LIRGs), a lower mean equivalent width 
(EQW6.2/im = 0.30 /im versus EQW6.2/^m = 0.55 /im). 
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and deeper average silicate absorption (sg.y^m = -1.05 
versus = -0.25). 

3) There is a general trend among the U/LIRGs for both 
silicate depth and MIR slope to increase with increasing 
L/jj. As L/fl increases, the temperature may rise as a 
consequence of the nuclei becoming more obscured and 
compact. As the dust temperature increases, the ris- 
ing portion of the blackbody emission spectrum shifts to 
shorter wavelengths, and warmer sources have increas- 
ingly more flux at 30 /xm, and thus steeper MIR slopes. 
The sources that depart from these correlations, in both 
cases, have very low PAH equivalent width (EQW6.2pm 
< 0.1 fim) consistent with their MIR emission being dom- 
inated by an AGN. 

4) Although less numerous (only 18% of the sample), 
LIRGs with the largest contributions from AGN (those 
with EQW6.2;im < 0.27 /im) cover a wider range of MIR 
slopes and silicate strengths than those sources of higher 
equivalent width that make up the majority of the sam- 
ple. The sources with extremely low PAH equivalent 
widths (EQW6.2pm<0.1 /xm) separate into two distinct 
types: relatively unobscured sources with a very hot dust 
component (and thus very shallow MIR slopes) and heav- 
ily dust obscured nuclei with a steep temperature gradi- 
ent. For the AGN-dominated LIRGs with low apparent 
obscuration, an upper limit to the MIR slope can be set 
at Fi^[30/im]/F^[15/xm]~4. The most obscured nuclei, 
however, have steeper MIR slopes due to most of their 
warm dust emission being hidden behind a large amount 
of cooler dust. 

suggesting 5) The LIRGs most likely harboring buried 
AGN (the obscured nuclei with S9.7^m<-1.75) all have 
EQW6.2/im <0.2 /im and lack any extended component 
to their MIR emission at 24 /im. Extreme levels of dust 
obscuration may simply be blocking PAH emission, or 
the higher densities in these nuclei may create an environ- 
ment where PAH dust grains are not present or the con- 
ditions are not appropriate for exciting them (i.e. lacking 
in photodissociation regions) . Their compact nature sug- 
gests that their obscuring (cool) dust is associated with 
the outer regions of the starburst and not simply a mea- 



sure of the dust along the line of sight through a large, 
dusty disk. 

6) U/LIRGs in the late to final stages of a merger have, 
on average, steeper MIR slopes and higher levels of dust 
obscuration. As merging galaxies coalesce and gas & dust 
is funneled towards the center, the nuclei become more 
compact and more obscured. As a result, the dust tem- 
perature increases leading also to a steeper MIR slope. A 
marked decline is seen for the fraction of high EQW (star 
formation dominated) sources as the merger progresses. 
The decline is accompanied by an increase in the fraction 
of composite sources while the fraction of sources where 
an AGN dominates the MIR emission remains low. 

7) Despite their dusty and starbursty nature, the aver- 
age nearby LIRG spectrum does not resemble the average 
composite (starburst + AGN) MIR spectrum from sub- 
millimeter galaxies at z^2. Both the average LIRG and 
ULIRG spectra are more absorbed at 9.7 /im and the av- 
erage LIRG has more PAH emission. However, once the 
AGN contributions are removed from the average LIRG 
and from the average SMG spectra, the PAH emission 
and level of silicate absorption of the average spectrum 
for starburst-dominate d SMGs (i.e. those without A GN 
spectral signatures; Me nendez-Delmestre et al.|2009 ) are 
fit well by the average starburst-dominated local LIRG. 

The Spitzer Space Telescope is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under NASA contract 1407. This research has made 
use of the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 
This research has made use of the NASA/IPAC Infrared 
Science Archive, which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics Space 
Administration. We would like to thank M. Cluver for 
many helpful discussions and K. Menendez-Delmestre 
for sharing her average SMG spectra. 
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